Does adsorption in a single nanogroove exhibit hysteresis? 
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A simple fluid, in a microscopic capillary capped at one end, is studied by means of fundamental 
measure density functional. The model represents a single, infinitely long nanogroove with long- 
range wall- fluid attractive (dispersion) forces. It is shown that the presence or absence of hysteresis 
in adsorption isotherms is determined by wetting properties of the wall as follows: Above wetting 
temperature, T w , appropriate to a single wall of the groove, the adsorption is a continuous process 
corresponding to a rise of a meniscus from the capped to the open end of the groove. For a sufficiently 
deep capillary the meniscus rise is shown to be a steep, yet continuous process taking place near 
the capillary condensation of a corresponding slit. However, for temperatures lower than T w the 
condensation exhibits a first-order transition accompanied by hysteresis of the adsorption isotherm. 
Finally, it is shown that hysteresis may occur even for T > T w as a consequence of prewetting on 
the side and bottom walls of the groove. 



The interfacial properties of fluids in confining geome- 
tries play a key role in several branches of physics and 
chemistry, and are vital for numerous engineering ap- 
plications, see, e.g., Ref [1 . Two of the most funda- 
mental interfacial phenomena involve a development of a 
liquid phase on a single wall and between two parallel un- 
bounded walls at a pressure below the phase coexistence. 
The first, known as complete wetting, is characterized by 
a growth of an adsorbed liquid film t ~ Ap~^° . Here i 
is the film thickness and A/i is a deviation of the chem- 
ical potential from its saturation value, /i sat (T), above 
a wetting temperature T w (which corresponds to a zero 
contact angle). The exponent (3^° is non-universal and 
depends on the character of molecular interaction; in par- 
ticular, f3g° = for short-ranged forces, while = 1/3 
when dispersion forces are involved [2 j. The second phe- 
nomenon, capillary condensation, differs from the previ- 
ous one in two respects. First, it is a first-order tran- 
sition, and, second, it reflects the finite size shift of the 
(3D) bulk phase boundary rather than surface phenom- 
ena at the wall. For a non-retarded van der Waals fluid- 
wall interaction and in the limit of a large distance L 
between the walls, the modified Kelvin equation [3l H] 
predicts the undersaturation at which the capillary con- 
densation occurs at T > T w : 



A/i(L) = /i sat - fJL C c(L) 
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(pi-Pg)(L-M) 



(1) 



Here 7 is the surface tension of a free liquid-gas interface, 
pi and p g are densities of the coexisting bulk phases, and 
£ is an adsorbed-film thickness. It should be noted that 
the details of the intermolecular forces in ([!]) are reflected 
only in the Derjaguin's correction (the factor of three 
multiplying the film thickness). 

While adsorption phenomena in systems that can be 
described by a one-dimensional density distribution are 
rather well understood, interfacial phenomena on pat- 
terned surfaces, which are characterised by a variation 



in the density distribution at least in 2D, have recently 
become the focus of considerable interest, see, e.g., Ref. 
[5]. The presence of a complex distribution of possible 
substrate topographies introduces a challenging task of 
relating geometrical and well-defined substrate models to 
the thermodynamics of fluid adsorption. Models suited 
to this task include linear wedges, cones, grooves, and 
pitted surfaces to name a few of the most popular [6]- 
[T3] . The interest in this research has been largely mo- 
tivated by new advances in nanofabrication techniques 
[T4] which made possible a direct comparison between 
theoretical predictions and experiment [15 . Apart from 
the theoretical interest, the adsorption in nanopatterned 
surfaces has attracted considerable attention due to its 
application in micro- and nanofluidics [T6] . 

The purpose of the paper is to give a microscopic de- 
scription of the adsorption of a simple fluid in a single, 
infinitely long groove of finite depth. In particular, we 
wish to know whether the adsorption isotherms in this 
model exhibit hysteresis or not. The main conclusion 
based on earlier studies [TTHT31 [T7| [T8] is that in the case 
of complete wetting the hysteresis vanishes due to the 
presence of the bottom wall. Here we present a detailed 
microscopic study of this issue and discuss the effect of 
the wetting properties on the character of the groove- 
adsorption. 

Consider a semi-infinite solid slab with a uniform one- 
body density p w spanning a domain S = R ® R (g) 
(— 00, L z ), L z > 0, in 3D Cartesians. Imagine that an in- 
finitely long groove of width L x and depth L z , occupying 
a subspace G = (0, L x ) <S> (—00, 00) (g) (0, L z ), is sculpted 
into the slab. We assume that the groove is subject to a 
potential: 



V(x,z) = 



x < o~ qn or x > L T — a in or z < <j ?/ 



V(x,z) , elsewhere , 



(2) 
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where V(x, z) — p w J w <fr w (\r — r'Ddr', W = S \ G, and 

Mr) = ~te w (^) 6 • (3) 

Thus, we assume that the fluid particles interact with 
the substrate atoms via long-ranged (dispersion) forces 
characterized by the parameters e w and <j w . 

After integration, V(x, z) can be expressed as follows: 

V{x, z) = Vxiz) + V 2 {x, z) + V 2 (L X -x,z), (4) 

with ViO) = ^ and V 2 (x,z) 

a w [i/j(x, z) + i/j(x, L z — z)) , where we have defined 

a™ = — \ire w p w (T% and 

,, N 2z 4 + xV + 2x 4 1 
^ z) = . - . 

The fluid-fluid interaction is given by 

{oo, r < a; 

-4e(^) 6 , a<r<r c ; (5) 
0, r > r c . 

In the following, parameters a and e will be used as 
length and energy units, respectively, and the cutoff is 
set to r c = 2.5 a. The wall parameters are fixed to 
<j w = <j and e w = 1.2 e. The wetting temperature for 
such a model is T w = 0.83 T c , where ksTc/e = 1.41 is 
the critical temperature. 

Within a classical density functional theory [19], the 
equilibrium density profile is found by minimizing the 
grand potential functional 

{l[p]=F[p} + Jdrp(r)[V(r)-p], (6) 

where \i is the chemical potential. J~[p] is the intrin- 
sic free energy functional of the one-body density, p(r), 
whose excess (over ideal gas) part is treated as a pertur- 
bation about a hard- sphere reference fluid: 

FeM = Fhslp] + \j drp(v) J dvp(v')4>(\v - r'|) . (7) 

The hard sphere part of the excess free energy functional 
is approximated by means of the Rosenfeld fundamental 
measure theory (FMT) [20] 

J=haW = i J dr^({n a }), (8) 

where $ is a function of the weighted densities {n a (r)}, 
and P = l/fc^T. Numerical calculations were carried out 
on a 2D grid with mesh 0.05 a. 

One should note that the functional ^ with JhsM 
obtained from the FMT is known to properly account 
for the short-ranged correlations that play a significant 
role in confined systems. In our nanogroove model one 



expects strong packing effects at the vicinity of the walls 
(as in the case of common slit models) but additional and 
particularly strong inhomogeneities are also expected at 
the edges of the groove. Finally, our DFT is thermody- 
namically consistent and yields the correct divergence of 
wetting film thickness for systems that exhibit dispersion 
forces. All these properties are of key importance for a 
reliable description of the phase behaviour of the model. 
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FIG. 1: Meniscus rise in a groove of width L x = 7 a (empty 
symbols) and L x — 12 a (full symbols) as a function of 
undersaturation (pb is the bulk density, p c is the critical 
density). The depths of the grooves are from the bottom 
L z /a — 20,30,40, and 50. The curves are the guides to the 
eye. The vertical lines denote densities corresponding to cap- 
illary condensation in an infinite slit of a width L x — 7 a 
(dashed) and L x = 12a (dotted). For T = 0.92 T c . 

In the following, we summarize the most-important 
features of adsorption in the model groove that can be in- 
ferred from our DFT calculations. We start by consider- 
ing temperature T = 0.92 T c , which is above the wetting 
temperature, T w , of a corresponding planar wall. In such 
a case, a meniscus separating a gas-like and a liquid-like 
phase is formed near the bottom end, and continuously 
rises upon increasing the chemical potential towards the 
two-phase coexistence region. Position of the meniscus 
can be defined as i(x) = J dz(p(x,z) — pb)/Ap, where 
Ap = pi — p g is the difference between liquid and vapour 
densities at the bulk coexistence. In Fig. 1 the menis- 
cus rise is shown for grooves of L x = 12 a and L x = 7 a 
and several depths. As is apparent, the groove filling is a 
continuous process, which should be contrasted to capil- 
lary condensation in infinite slits. Furthermore, as L z is 
increased, the meniscus rise becomes steeper and steeper 
upon approaching the capillary condensation indicated 
by the vertical lines. Such results should be compared 
with a schematic plot (Fig. 1) in reference [T3] , 

Next consider T = 0.81 T c , which is slightly below 
T w . Now, the scenario of the fluid adsorption dra- 
matically changes. As displayed in Fig. 2, the sys- 
tem exhibits a hysteresis in adsorption per unit length 
r = J G p{x, z) dx dy, indicating a first order transition 
resembling capillary condensation in infinite slits. In or- 
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FIG. 2: Adsorption isotherm at temperature T — 0.81 T c in 
a groove of width L x = 12 a and depth L z = 20 a; pb is the 
bulk density, p g is the density of a saturated gas. Full sym- 
bols denote the adsorption path, empty symbols denote the 
desorption path. The vertical dashed line denotes a density 
where two distinct states coexist. 



der to understand the different behaviour of the adsorp- 
tion in the two cases, let us discuss the mechanism of the 
groove filling for T > T w and T < T w . 
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FIG. 3: Schematic illustration of the slab model. 

For T > T w , the groove filling can be described with 
the help of a simple slab model, as schematized in Fig. 3. 
Substituting a sharp-kink approximation for the density 
profile in (|6| , the excess (over the groove completely filled 
with a liquid) grand potential: 



L 

+7[2(iz - 



- = {vt -P){Lx ~ 
4) + (L x - 24)] 



24) (L x -4) (9) 
+ Pjl^PL f v(r)dr , 



where pf is the bulk pressure of the metastable liquid 
and p is that of the gas reservoir. When minimized with 
respect to £ z (minimization wrt to £ x reproduces Eq. ([!])), 
one obtains in the limit of large L z : 



Ap- 



A 



■0(1/0 = 0, 



where Ap 



Pi 
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9 a w (pi—p g )L x (L x — 2£ x ) 
8 L x -3£ x 



L x — 3£ x 

Comparing diok 



(10) 

and A = 
with (!]) 



and assuming that Afi(L x ) is small, the meniscus rise 
satisfies (cf. Ref. [T3] ) 

-1/4 



4 ~ (Vcc(L x ) - 



(11) 



The same conclusion was obtained in Ref. [13]. However, 
in Ref. [13] the authors ignored the influence of £ x , while 
the connection with Eq. has been made more explicit 
here. Thus, for T > T w one observes a continuous filling 
of a groove, similar to a complete wetting on a planar 
wall. However, the differences with the latter are: i) the 
deep groove becomes filled in the limit \i — >• /jl(L x ) rather 
than for \i — >• /i sa t5 an d ii) the critical exponent of the 
groove filling is 1/4, i.e., the process is somewhat slower 
than the complete wetting. 

In contrast, for T < T w only a microscopically thin 
film is formed on a free wall. Therefore, the system be- 
haves in a qualitatively same fashion as in an open slit 
and undergoes a first-order transition accompanied by 



pronounced hysteresis. We emphasize that Eq. 10) can- 
not be applied anymore in this regime, since inclusion of 
higher order terms due to repulsive interactions would be 
needed. These would give a global minimum of ^ ex at a 
finite distance of £ z . 

Finally, let us consider temperature T = 0.91 T c . Such 
a temperature is again above T w but this time (slightly) 
lower than T cs , the prewetting critical temperature cor- 
responding to a single wall. It is well known [21] that the 
prewetting transition between a thin and thick layers on a 
single wall may also occur in slits, but the phenomenon is 
preceded by capillary condensation unless the separation 
between the walls is large. It is interesting to examine 
how prewetting is reflected in our model groove in the 
absence of capillary condensation (above T w ). In Fig. 4 
an adsorption isotherm for a model L x = L z = 50<j is 
depicted. We observe that the otherwise-continuous in- 
crease in adsorption, due to meniscus rise, exhibits two 
small van der Walls loops, both slightly below the bulk 
density appropriate to prewetting. The first of the two 
loops corresponds to a discontinuous jump in the menis- 
cus height, reflecting prewetting on the bottom wall, cf. 
Fig. 5 (top). Such a transition is shifted to a higher un- 
dersaturation compared to prewetting on a single wall, 
pjTVpsat, ^ ue t° the influence of the side walls that de- 
creases the potential inside the groove. The second loop 
corresponds to prewetting on the side walls, cf. Fig. 5 
(bottom). Again, this transition is shifted below p% w /p s &t 
but the difference is smaller, since the effect of the bot- 
tom wall on the vertical liquid films is relatively weak (for 
large z the potential of the bottom wall decays as z~ 3 ). 
However, the scenario may be different and the sequence 
of the two transitions be reversed if the side walls are 
stronger adsorbents than the bottom wall. 

In summary, we have employed a mean-field FMT- 
DFT to describe adsorption in a single nanogroove with 
attractive walls that exhibit a first-order wetting transi- 
tion. At high temperatures, the results support the pic- 
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FIG. 4: Adsorption isotherm at temperature T — 0.91 T c in 
a groove of width L x = 50a and depth L z — 50 a. Black 
symbols denote the adsorption path, red symbols denote the 
desorption path. The lines are the guides to the eye. The 
vertical dashed line denotes a density corresponding to the 
prewetting transition in an open slit of width L = 50cr. 
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FIG. 5: Density distributions (x-z projection) in a groove with 
attractive walls at temperature T — 0.91 T c . L x — L z — 50a. 
Top: coexisting states at a density pb/p g = 0.93695. Bottom: 
coexisting states at a density pb/pg — 0.94813. 



ture of the onset of a liquid phase due to a heterogeneous 
nucleation on a bottom wall, followed by a continuous 
rise of the meniscus separating vapour-like and liquid-like 
phases, and no hysteresis in adsorption appears in this 
case (see Fig 1). However, as we have demonstrated, such 
a scenario is valid only for temperatures above the wet- 
ting temperature of a corresponding single wall. Below 
T W1 the adsorption exhibits a first-order transition (ac- 
companied by hysteresis) as in an open slit. Thus, T w , 
which represents a crossover between partial and com- 
plete wetting for a single wall, represents a boundary 
separating first-order and continuous transition regimes 
for a nanogroove. Some care is needed, though. For 
large widths, a wedge filling may take place at the cor- 
ners of the groove at temperatures above a filling tem- 



perature, which may cause a shift of this boundary below 
T w ; for the model considered here the filling temperature 
(of a rectangular wedge) is Tf « 0.76 T c . We haven't 
observed such a phenomenon even for a groove width as 
large as L x = 50 a, but one cannot rule out this possibil- 
ity. Finally, it has been shown that in the temperature 
interval (T W ,T SC ), the groove-adsorption isotherm may 
exhibit two loops, reminiscent of prewetting. The first 
(if the bottom and side walls are adsorbents of the same 
strength) corresponds to a jump of the meniscus, which 
is followed by an discontinuous thickening of the liquid 
films on the side walls. 

We believe that the observed phenomena will moti- 
vate further experimental measurements and give some 
contribution into the long-standing issue of a well-known 
but still actively debated phenomenon of hysteretic be- 
haviour of fluid adsorption in mesoporous media. 
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